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The dynamics of vortex lattices in stirred Bose-Einstein condensates have been studied at finite 
temperatures. The decay of the vortex lattice was observed non-destructively by monitoring the 
centrifugal distortions of the rotating condensate. The formation of the vortex lattice could be 
deduced from the increasing contrast of the vortex cores observed in ballistic expansion. In contrast 
to the decay, the formation of the vortex lattice is insensitive to temperature change. 
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Gaseous Bose-Einstein condensates (BEC) have be- 
come a testbed for many-body theory. The properties 
of a condensate at zero temperature are accurately de- 
scribed by a nonlinear Schrodinger equation. More re- 
cently, theoretical work on the ground state properties 
of condensates Q has been extended to rotating con- 
densates containing one or several vortices and their 
dynamics. These include vortex nucleation 
tallization of the vortex lattice ||, and decay 
perimental study has focused mainly on the nucleation 
of vortices, either by stirring condensates directly with a 
rotating anisotropy |^-[Io[| or creating condensates out of 
a rotating thermal cloud Here we report the first 

quantitative investigation of vortex dynamics at finite 
temperature. The crystallization and decay of a vor- 
tex lattice have been studied and a striking difference is 
found between the two processes: while the crystalliza- 
tion is essentially temperature independent, the decay 
rate increases dramatically with temperature. 

The method used to generate vortices has been out- 
lined in previous work P,^2[. Condensates of up to 
75 million sodium atoms (> 80% condensate fraction) 
were prepared in a cigar-shaped Ioffe-Pritchard mag- 
netic trap using evaporative cooling. The radial and 
axial trap frequencies of u> x = 2ir x (88.8 ± 1.4) Hz, 
Lu y = 2ir x (83.3 ±0.8) Hz, and u z = 2ir x (21.1 ±0.5) 
Hz, respectively, corresponded to a radial trap asymme- 
try of e r = (lu 2 x - u>l)/{u)l + uj 2 y ) = (6.4 ± 0.2)% and an 
aspect ratio of A = u) x /ui z — (4.20±0.04). The relatively 
large value of the radial trap asymmetry is due to grav- 
itational sag and the use of highly elongated pinch coils 
(both estimated to contribute equally). The radio fre- 
quency used for evaporation was held at its final value to 
keep the temperature of the condensate roughly constant 
throughout the experiment. The condensate's Thomas- 
Fermi radii, chemical potential, and peak density were 
R r = 28 fjtm, R z = 115 fim, 300 nK, and 4 x 10 14 cm~ 3 , 
respectively, corresponding to a healing length £ ~ 0.2 
/xm. 

Vortices were produced by spinning the condensate for 
200 ms along its long axis with a scanning, blue-detuned 
laser beam (532 nm) [pUlq] . For this experiment two sym- 



metric stirring beams were used (Gaussian waist w = 5.3 
/xm, stirring radius 24 fim). The laser power of 0.16 mW 
per beam corresponded to an optical dipole potential of 
310 nK. After the stirring beams were switched off, the 
rotating condensate was left to equilibrate in the static 
magnetic trap for various hold times. As in our previous 
work, the vortex cores were observed using optical pump- 
ing to the F = 2 state and resonant absorption imaging 
on the F = 2 to F — 3 cycling transition After 42.5 
ms of ballistic expansion the cores were magnified to 20 
times their size, £, in the trap. 

The decay of the vortex lattice can be observed by al- 
lowing the condensate to spin down in the trap for vari- 
able times before the ballistic expansion and subsequent 
observation of the vortex cores. However, this method 
requires destructive imaging and suffers from shot-to- 
shot variations. In situ phase-contrast imaging is non- 
destructive, but the vortex cores are too small to be re- 
solved. Instead, we monitored the centrifugal distortion 
of the cloud due to the presence of vorticity. This is a 
quantitative measure for the rotation frequency, Q, of the 
lattice, and therefore the number of vortices Such 
distorted shapes have been observed previously for rotat- 
ing clouds |,0Jl|. 

The shape of a rotating condensate is determined by 
the magnetic trapping potential and the centrifugal po- 
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where M is the atomic mass. From 



the effective radial trapping frequency u> r — ► uj' r = 
■ f2 2 , one obtains the aspect ratio of the rotating 



y/tL 

cloud, A' — uj' r /u) z , 



A' = Ay/1 - {VL/LO r f 



(1) 



For the decay measurements, the condensate was 
stirred for 200 ms, producing ~ 130 vortices on aver- 
age at the coldest temperature (determined using time- 
of- flight imaging). After the drive was stopped the cloud 
equilibrated in the stationary magnetic trap. Ten in situ 
images of each condensate were taken at equal time in- 
tervals using non-destructive phase contrast imaging de- 
tuned 1.7 GHz from resonance (see Fig. [l]a). We verified 
that the decay was not affected due repeated imaging 
by varying the time between exposures and observing no 
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difference. The theoreticai treatment of the decay has 
two limiting cases, one where the thermal cloud is nearly 
stationary the other where the thermal cloud is closely 
following the rotating condensate Because the ther- 
mal cloud could not easily be discerned in images taken 
with 1.7 GHz detuning, we used light closer to resonance 
(400 MHz detuning) to determine the shape of the ther- 
mal cloud (Fig. |l]b). 
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fits to time-of-flight images of the condensate before the 
stirring. The temperature was derived from these values 
using the scaling theory for an interacting Bose gas [jl5| . 
At a fixed rf frequency for evaporation, the centrifugal 
potential lowers the trap depth by a factor 1 — (il/ui r ) 2 
for a thermal cloud rotating at frequency £1. Evapora- 
tion should lower the temperature of the rotating cloud 
by the same factor. 

Fig. ||a shows how the aspect ratio of the condensate 
approaches its static value as the vortex lattice decays. 
The decrease of the angular speed appears to be expo- 
nential (Fig. ||b) and strongly depends on temperature 
(Fig. 1). 




FIG. 1. a) Spin down of a rotating condensate in a static 
magnetic trap. The first phase contrast image is taken 100 
ms after turning off the drive, with each subsequent image 
spaced by 100 ms. The rotating vortex lattice caused a radial 
swelling of the condensate, reducing the aspect ratio. As vor- 
tices leave the system the aspect ratio approaches its static 
value of 4.20. The field of view is 25 ^tm by 75 fim. b) Ob- 
servation of a rotating thermal cloud. The parameters are 
identical to (a), but the probe light detuning is closer to res- 
onance, enhancing the sensitivity to the more dilute thermal 
cloud. The phase shift of the dense condensate exceeds 2n 
and is displayed as saturated in the image. The apparent loss 
in atom number is due to Rayleigh scattering of the probe 
light, c) The inner and outer contours represent the aspect 
ratio of the condensate and thermal cloud, respectively, as 
obtained from two-dimensional fits to phase-contrast images. 

The aspect ratio of the condensate was obtained by 
fitting a 2D Thomas-Fermi profile to the 1.7 GHz images. 
The aspect ratio of the thermal cloud was determined 
from the 400 MHz images by masking off the condensate 
in the inner region of the cloud, and then fitting a 2D 
Gaussian to the image. The comparison of the contours 
of the condensate and thermal cloud (condensate fraction 
of N c /N — 0.62) shows that the thermal cloud is also 
rotating. The images for the thermal cloud show that 
the aspect ratio decays from 3.0 to 4.2 on the same time 
scale as the condensate. This corresponds to an initial 
rotation rate of ~ 2/3 that of the condensate. 

The damping rate of the vortex lattice, obtained from 
exponential fits to the data, was studied at different tem- 
peratures by varying the condensed fraction of atoms in 
the trap. The condensate fractions were obtained from 
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FIG. 2. Decay of a vortex lattice at finite temperatures, 
a) The aspect ratio of the condensate approaching its static 
value. Each point represents the average of 10 measurements. 
Error bars given for (•) are statistical and are typical of all 
four data sets, b) The decay of the rotation rate of the lattice 
for the same data using Eq. (|l|). The decay of the rotation 
depends strongly on the thermal component. 

In addition to the decay process, the formation of the 
vortex lattice has also been examined. After rotating the 
condensate it typically took hundreds of ms for the lat- 
tice to form (see Fig. 4 in ref. ||). One may expect 
the lattice to already form in the rotating frame dur- 
ing the stirring because the lattice is the lowest energy 
state for a given angular momentum. This absence of 
equilibration in the rotating frame is presumably due to 
heating and excitation of collective lattice modes by the 
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stirring beams. The crystallization of the vortex lattice 
was studied by determining the contrast or visibility of 
the vortex cores as a function of equilibration time. To 
avoid any bias, we used an automated vortex recognition 
algorithm. First, each image was normalized by divid- 
ing it by a blurred duplicate of itself. A binary image 
was then created with a threshold set to the value at 
which the algorithm reliably detected almost all vortices 
that were identified by visual inspection of equilibrated 
images (there was less than 5% discrepancy). Clusters 
of contiguous bright pixels within a circular area were 
counted as "visible" vortices. Fig. [|a shows three vortex 
lattices after different equilibration times with the visible 
vortices identified. 
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FIG. 3. Decay rates for vortex lattices at several tempera- 
tures. The rates are determined from data taken during the 
1st second of equilibration. Each point represents the average 
of 5 measurements. The inset figure emphasizes the strong de- 
pendence of the decay rate on temperature, T. The transition 
temperature was T c = 790 nK and varied by less than 7 % 
for the data shown. 

The lowest temperature (337 nK, corresponding to a 
condensate fraction of w 80%) was close to the chemical 
potential. Further cooling would have resulted in smaller 
condensates. The highest temperature (442 nK) was lim- 
ited by the rapid decay of the vortex lattice. 

Fig. ^b,c show that the formation of the lattice depends 
very weakly on temperature, if at all, over the range stud- 
ied. The larger number of vortices crystallized at lower 
temperatures (Fig. ^b) is due to the strong temperature 
dependence of the lattice decay rates, which differed by 
more than a factor of 5. We can correct for this by esti- 
mating the number of vortices N v (t) as a function of time 
from the centrifugal distortion measurements described 
above as N v (t) = 2K,fl(t)7rR(t) 2 , where n = h/M is the 
quantum of circulation, M the sodium mass and R the 
radial Thomas-Fermi radius. By normalizing the number 
of visible vortices by this estimate, we deduce the vortex 
visibility as a function of time. These lattice formation 
curves overlap almost perfectly for different temperatures 
(Fig.|c). 
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FIG. 4. Crystallization of the vortex lattice, a) The top 
row shows three condensates that have equilibrated for 50, 
150 and 300 ms, respectively, and have 48, 86 and 140 vor- 
tices recognized as "visible" by our algorithm. The bottom 
row shows the same condensates with the "visible" vortices 
circled. The field of view was 1.4 mm by 1.6 mm. b) Growth 
of the number of visible vortices for several temperatures ex- 
pressed by the condensate fraction (N c /N). c) Visibility of 
vortices derived from the data in figure b, normalized by the 
number of vortices inferred from centrifugal distortion mea- 
surements. 



The decay of a vortex lattice is discussed in refs. 
The process is modelled as a two-step transfer of angular 
momentum. The rotation of the condensate is damped 
due to friction with the thermal cloud (fi c -th), which 
is in turn damped via friction with the trap anisotropy 
(fith-c,.)- Even at our highest temperatures, fi c -th is suf- 
ficiently larger than [ith-e r to spin up the thermal cloud 
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to ~ 2/3 of the lattice rotation rate. fJ,th-e r is calculated 
in ref. |ll| for a classical Boltzmann gas. For our pa- 
rameters the relaxation time r rot for the thermal cloud 
rotation is r rot = 4r, where r is the relaxation time for 
quadrupolar deformations, r is approximately 5/4 of the 
time between elastic collisions r e i = 2/nv t h& [0, where 
n is the atomic density, Vth the thermal velocity and a 
the elastic cross section. For a quantum-saturated ther- 
mal cloud, this gives a relaxation time proportional to 

The presence of the condensate adds additional mo- 
ment of inertia, extending the rotational relaxation time 
by a factor / in ortia = (ith + Ic)/hh, where I th and I c 
are the moments of inertia of the thermal cloud and the 
condensate, respectively. At high temperatures this fac- 
tor approaches one, because the moment of inertia of the 
thermal cloud is dominant. For very low temperatures 
/inertia ~ (i c I Ith)- This ratio is §] 



Ic_ = 15C(3) N c (i 
J th 16C(4) N th k B T 



(2) 



where the numerical pre- factor is 1.04. Thus, for a non- 
interacting Bose gas in the low temperature limit /i ne rtia 
scales as T~ A and the relaxation time of the vortex lattice 
should scale as T~ 6 . The observed temperature depen- 
dence of the relaxation time, decreasing by a factor of 17 
for a 60% increase in the temperature, agrees fortuitously 
well given the approximations made in the theory. The 
absolute relaxation rates predicted by refs. |j 16 are on 
the order of 1 s _1 , in reasonable agreement with our re- 
sults. However, for a quantitative comparison to theory 
it will be necessary to further characterize the rotating 
thermal cloud, which may not be in full equilibrium. 

In contrast to the decay, the crystallization process of 
the vortex lattice was essentially independent of temper- 
ature. This was unexpected because all dissipative pro- 
cesses observed thus far in BECs, including the decay of 
the vortex lattice, have shown a strong temperature de- 
pendence [tL8|JlS|] . D. Feder has numerically simulated a 
stirred condensate using only the Gross-Pitaevskii equa- 
tion. The results show that the absence of dissipation 
leads to rapid, irregular motion of the vortices. The ad- 
dition of a phenomenological dissipative term in ref. j^] 
resulted in the formation of triangular vortex lattices. 
However, the origin of the dissipation was not identified. 
Our results suggest that the time-limiting step for the 
evolution of a vortex tangle into a regular lattice does 
not strongly depend on temperature. One possibility is 
that the thermal cloud is not directly involved. This may 
be similar to the reconnection of vortices or the damp- 
ing of Kelvin modes, where the spontaneous emission of 
phonons act as a dissipative mechanism |^,^J . Another 
possibility is that the rearrangement of the vortices into 
a rectilinear lattice is slow and not limited by dissipation. 

In conclusion, we have studied the crystallization and 
decay of vortex lattices. Both processes are dissipative 



and require physics beyond the Gross-Pitaevskii equa- 
tion. The dynamics of vortices nicely illustrates the in- 
terplay of experiment and theory in the field of BEC. 
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